Convective instabilities in two superposed horizontal liquid layers heated laterally 
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' This work is devoted to the theoretical study of the stabihty of two superposed horizontal liquid 

, layers bounded by two solid planes and subjected to a horizontal temperature gradient. The liquids 

are supposed to be immiscible with a nondeformable interface. The forces acting on the system 
{^jq' are buoyancy and interfacial tension. Four different flow patterns and temperature profiles are 

r-{ ^ found for the basic state. A linear perturbative analysis with respect to two and three dimensional 

perturbations reveals the existence of three kind of patterns. Depending on the relative height of 
both liquids several situations are predicted: either wave propagation from cold to the hot regions, 
or waves propagating in the opposite direction or still stationary longitudinal rolls. The behavior 
CNj , of three different pairs of liquids which have been used in experiments on bilayers under vertical 

gradient by other authors have been examined. The instability mechanisms are discussed and a 
qualitative interpretation of the different behaviors exhibited by the system is provided. In some 
configurations it is possible to find a codimension-two point created by the interaction of two Hopf 
modes with different frequencies and wavenumbers. These results suggest to consider two liquid 
layers as an interesting prototype for the study of propagation and interaction of waves in the 
^ , context of the Benard-Marangoni problem. 
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I. INTRODUCTION 



There has been recently an increasing interest in the study of buoyant-thermocapillary flows induced by horizontal 
thermal gradients [J. Such situations occur in oceanographic currents due to non homogeneous heating, and their 
interplay with the atmosphere gives rise to complex meteorological phenomena. They are also important in crystal 
QQ ■ growth. The floating zone process generates thermocapillary convective cells in the region of crystallization which 

\ require better understanding and control. 
00 ■ When a horizontal temperature gradient is acting on a horizontal liquid layer enclosed between two horizontal walls 
' the problem is much more complicated than in the the more classical case of vertical heating. First, the reference state 
, is no longer at rest, giving rise to a general flow and a non-linear vertical temperature profile. A non trivial problem 
will therefore consist in determining the basic reference temperature and velocity fields. Another peculiarity is that 



^ the threshold of instability depends on the Prandtl number Pr. The problem of thermocapillary instabilities induced 
, ^ by a lateral heating was first studied by Smith and Davis ■ They predicted the presence of hydrothermal waves and 
^ stationary rolls. More details are provided in a subsequent work by Smith Isj , who showed the existence of two different 



mechanisms of instability depending on Pr. At low Pr, the energy necessary to sustain the disturbances comes from 
the horizontal applied temperature field and the hydrothermal waves propagate in a direction perpendicular to the 
horizontal temperature gradient. At high Pr, the energy is extracted from the vertical temperature field by vertical 
>^ ' convection and it gives rise to hydrothermal waves propagating parallel to the temperature gradient. At intermediate 
, Pr, the mechanism is a combination of the previous effects and leads to waves forming an angle with the streamwise 
\^ • direction. Parmentier et al. Q studied the coupling of buoyant and thermocapillary driven instabilities, in the 
case of systems with lower and upper insulating bounding surfaces within a 3D linear formulation. Their numerical 
results display the three kinds of behavior described above. For thermally conducting surfaces, Gershuni et al. Q 
find stationary rolls for Pr > 1. Other theoretical works have mainly focused on comparison with experiments. 
Mercier and Normand |^ performed a linear analysis to explain the experimental results by Daviaud and Vince 
for a silicon oil with Pr = 10. They notice that the transition between traveling waves and stationary rolls 
observed when the liquid depth is increased depend on the heat transfer between the liquid and its environment. In 
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their discussion about absolute, global and convective instabilities Priede and Gerbeth show that the threshold of 
absolute instability fits better the experimental data of 1 cSt silicon oil {Pr = 13.9) |9| than the value predicted from 
convective instability. A non-linear analysis of stability for purely thermocapillary convection has been carried out by 
Smith |lOl| , who determines the range of existence of the two oblique traveling waves predicted by the linear theory. 
More recently, new results have been reported pj| on hydrothermal waves in different geometries, measuring their 
frequency, wavenumber and angle of propagation. 

In the previous theoretical works it is assumed that the liquid is surmounted by a gas phase above that remains 
passive (one-layer approximation). In some cases 0,1^] the effect of the adjacent gas is modelled by a Biot number, 
a phenomenological parameter that characterizes the heat transfer between liquid and gas. The Biot number was 
originally defined for steady conditions; however when convection takes place, the dynamics of thermal and mechanical 
perturbations in the gas may have a strong influence on the convection inside the liquid. In such a case, the introduction 
of a Biot number becomes questionable and the one-layer approximation should be revised. In the general case, a full 
two-layer hydrodynamical description is required. Furthermore well controlled experiments involve always two fluid 
layers 14]. Such systems are typical of many practical situations as Earth's mantle convection or encapsulated crystal 
growth. 

In the past, most theoretical [3| and experimental studies about Rayleigh-Benard convection in two 

immiscible liquids have been devoted to vertical heating. Only few contributions are found dealing with a horizontal 
temperature gradient. To our knowledge the first work on the subject is by Villers and Flatten they measured the 
velocity profiles in each layer as a function of the height in a system formed by water and heptanol. They also developed 
a simple theoretical model to calculate the horizontal velocity profile as a function of the ratio of viscosities, expansion 
coefficients and thicknesses of both layers. Later on, Doi and Koster j3| studied theoretically the thermocapillary 
convection under microgravity conditions in two immiscible liquid layers with a free upper surface, in 2-D geometries. 
They determine under which conditions the lower layer remains at rest. Moreover they carry out numerical simulations 
in a box of aspect ratio 4 in order to analyze the effect of the vertical walls. Numerical simulations in cavities of 
different aspect ratios for coupled thermocapillary and buoyancy-driven convection were performed by Liu et al. , 
and an asymptotic solution for the velocity in the limit of infinite aspect ratio has also been derived . 

The objective of the present work is to study coupled thermocapillary and buoyancy convection in a two-layer system 
of infinite horizontal extension subject to lateral heating and to find out the regimes occurring in this configuration. 
A linear approach to the problem is carried out, showing the kind of oscillatory and stationary behaviors that occur 
in this system. 

The paper is organized as follows: In section^ the general equations governing the problem are established. Section 
mil is devoted to the derivation of the temperature and velocity profiles of the basic fiow, the explicit expressions are 
reported in the Appendix. In section lTVl is carried out a linear stability analysis of the basic state towards bidimcnsional 
and three dimensional perturbations. Conclusions are drawn in Section ! VII 



II. PROBLEM FORMULATION 

We consider a system of two horizontal superimposed immiscible liquid layers of thickness h^^\ densities p^'^\ 
kinematic viscosities i>^-^\ coefficients of volume expansion a^-^\ thermal conductivities A^'^* and thermal diffusivities 
K*-*-*, superscript i — 1,2 refers to the lower and upper fluid respectively. The system is inflnite in the horizontal 
extension and is limited in the vertical direction by two horizontal rigid and thermally conductive walls (see Fig. 

The interface between both liquids is assumed to be horizontal and nondeformable. The origin of the Cartesian 
reference system is fixed at the interface. 

The system is heated, with a constant temperature gradient /3 imposed along the horizontal direction, which 
produces a conducting temperature profile given by T = r_ -I- (3x, where T_ is the temperature of the cold side. 
The Boussinesq approximation is taken for granted in both layers. The system is subject to the gravity field and its 
density is given by the state equation p^'^ — p^^ [1 — a^*^ (T — r„)]. The interfacial tension which acts at the interface 
admits a temperature variation approximated by the linear equation of state a = ao — 7(T — T_) where ^ — is 
a positive constant (at least for the pairs of liquids considered in this work). 

The governing equations expressing the balance of mass, momentum and energy in the Boussinesq approximation 
are given by 
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V-u« = 0, (la) 
-— + (u«-V)uW = --ll^-g[l-a(TW-r_)]z + i/«V2u«, (lb) 

-— + (u«-V)T« = K^V^rW, (Ic) 

wherein u^^^^ {u''^\ v^'^\ w^^^ ) is the velocity field, p^'-* the pressure, g the acceleration of the gravity, and z the unit 
vector in the vertical direction. The boundary conditions at the rigid and thermally conductive bottom and upper 
walls are 



z 

z 



= ^ u(i) = 0; T(i) = (T^ = r_ + Px), (2a) 

^ ^(2) ^ Q. y(2) ^ = r_ + /3a;). (2b) 



Furthermore, one must include the conditions of continuity of the temperature, heat flow and velocity at the interface 
(the normal velocity component is zero because the interface is supposed to be nondcformable) 

A(i)9,T(i) = A(2)9,T(^), a.u'^' - 9,u« = -^V,t('\ (3) 

the last relation in ^ expresses the balance between the tangential stresses at the interface. To transform the 
governing equations and boundary conditions in dimensionless form, the following scales are selected: for length h^^\ 

time h^^^ / velocity K'-^'>/h^^\ pressure p'-^^ v^^'' k'-^^ / h^^'' and temperature Ph'^^\ In terms of non-dimensional 
quantities, the balance Eqs. for the lower liquid read as 

V-u(i) = 0, (4a) 



Pr-'[^ + (u(i) . V)u(i)] = -Vp(i) - MT^'^ T^)]^ + V^u(i), (4b) 



dt 

while for the upper liquid 



+ (u(i) • V)r(i) = V^T^i), (4c) 



V-u(2) = 0, (5a) 

(5b) 



+ (u(2) . V)r(2) ^ Kv2r(2). (5c) 
The nondimcnsional boundary conditions are: 

z = -1 ^ u(i) = 0, T(i) = T„, (T^ = r_ + a;); (6a) 

z = a ^ u(2) = 0, r(2) = (T^ = T_ + x); (6b) 

Z = 0^ U(1)=U(2), = ^ 0, r«=r(2)^ (g^) 

In (|4I6|I the following nondimcnsional parameters have been introduced a — a^'^^ /a'^^\ n — kP'^ / v — u'^'^^ /v'^^\ 
A A'^V^'^'; a = p = ^ ^ p(2)^{2) I p{i)^(i) ^ rpj^g Prandtl number is defined with respect 

to the liquid 1, i.e. Pr = i/'-^^/k'^', as well as the Rayleigh number Ra = - — ^ and the Marangoni number 

Ma - --^ ^(1)^(1,^(1) • 
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III. THE BASIC STATE 

When a horizontal gradient is imposed, a stationary basic flow sets in each hquid with a horizontal velocity 
component depending on the vertical coordinate u*^*^ — (uq*' (z), 0, 0). On the other hand, the basic temperature 
profile is the superposition of the imposed horizontal gradient and the vertical profile t'^'^^{z) generated by the fluid 
motion: T^*) = T_ + a; + r'^*)(z) (in nondimensional form). To calculate u"^ [z) and r(*)(z) we replace 

uW = (mW(z),0,0), T« -r_ =a; + r(*)(z), (7) 

in I0J and and eliminating the pressure one obtains the following equations governing the behavior of the velocity 
and temperature fields in both liquids: 

5^3 = Ra , 9^2 t(i) = (8a) 

vdz^u"^^ — aRa , kOz^t^^^ ~ Uq^\ (8b) 

with the corresponding boundary conditions 

z = -1 ^ = 0' 'T^^^ = 0; (9a) 

z^a^ = 0, r(2) = 0; (9b) 

z = 0^ r(i)=T(2), 4')=4'\ (9c) 
9,t(i) = Aa,T(2), - dj^'^ = Ma. 

The return flow condition which requires that the net flow through a vertical plane must vanish in each layer, i.e. 

\«=o , r^'^^o, (10) 

-1 JQ 

provides two additional conditions that allows us to calculate explicitly the velocity and temperature fields Uq and 
T^^\ whose expressions in terms of z are found in the Appendix. 

As shown by Eqs. H14I17|I the basic velocity and temperature profiles are the sum of a thermocapillary term 
proportional to the Marangoni number Ma and a buoyancy term proportional to the Rayleigh number Ra. The 
velocity profile is a second order polynomial in z in the interfacial term and of third order in the buoyancy. Whereas 
the temperature profile is fourth order in the interfacial term and fifth in the buoyancy. 

The velocity and temperature profiles depend on the ratios of the various transport coefficients and depths. For a 
given experimental set-up the heating does not affect the shape of the basic velocity and temperature profiles since 
the temperature gradient appears as a constant factor in Uq^ and t'*^ through the Rayleigh and Marangoni numbers. 

Each velocity profile has a root located at its corresponding rigid wall (as imposed by the boundary conditions 
H6al6b|l ^. There is also a second root inside each layer due to the return fiow condition. In absence of gravity, this 
root is located at z = —1/3 in the lower layer and at z = a/3 in the upper fluid, i.e. one third of the depth counted 
from the interface, and the flow consists of one convective cell in each layer. When gravity is taken into account the 
position of these roots depend on the properties of the fluids and a third root can be present, allowing for a second 
convective cell in each layer. The temperature profiles have also a root located at the bounding horizontal walls. In 
presence of gravity it is possible to find until five roots in each layer. 



A. Velocity and temperature profiles 

The flow is driven by interfacial tension gradients, and by density differences when the gravity is acting. As usually, 
it is assumed that the surface tension decreases with temperature (7 > 0), so that a horizontal temperature gradient 
gives rise to tangential forces that drive the fluid from hot to cold regions. Therefore, with only thermocapillary 
effects present, a general circulation around the interface is established that will drag the fluid from the hot to the 
cold side. In a finite container, by continuity, the fluid in the upper layer raises near the cold side and falls down 
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near the hot one; on the contrary, in the lower layer it falls along the cold side and raises up at the hot one. One 
would therefore observe two counter-rotating cells. Buoyancy forces also drive the fluid from hot (lower density) to 
cold regions (higher density). Since the fluid raises along the hot side and falls at the colder one, in both cells, the 
buoyancy will favor the formation of two co-rotating cells. In coupled thermocapillary-buoyancy convection, buoyancy 
and interfacial forces are acting in the same direction in the lower layer, just like in one-liquid systems with 7 > 0. 
However they act in opposite directions in the upper layer, where the forces are competing (see Fig. I^J. Thus the 
liquid 2 exhibits a scenario similar to that of a one-liquid system with 7 < . (This occurs in some ceramics and 
liquid alloys, like Ag Pb, |2ll|. ) 



In Figs. 3(a) and 3(b)| are displayed the basic velocity and temperature proflles. Along the horizontal axis is 
reported the total depth (/ly) of the two superposed fluids while the vertical axis gives the percentage of depth of 
the lower liquid with respect to the total depth, i.e. the fractional bottom depth hi {hi = 100/(1 -|- h'''^^ /h'^^^)). The 
curves refer to a configuration formed by perfluorinated IIT-70 (lower fluid) and 5 cSt silicon oil (upper fluid), this 
pair of fluids has been used in a recent experiment with vertical heating by Juel et al. jl5| . (The physical properties of 
these liquids are listed in Tab. HI) In both fi gures, four different regions denoted (I) to (IV) are distinguished. Let us 



pay attention to Fig. 3(a) The region (I) extends to all hr and covers principally high values of hi. The basic state 
consists of a clockwise convective cell in liquid 2, and a counterclockwise cell in liquid 1; in this region the direction 
of the general circulation is governed by the interfacial forces. In region (II), the competition between buoyancy and 
interfacial forces gives birth to two sublayers in liquid 2. The motion in the sublayer close to the interface is dominated 
by the interfacial forces, but near the top wall it is the buoyancy that contributes to create a second sublayer. In 
regions (I) and (II), the general flow structure in liquid 1 is not changed by the motion of the upper liquid. However 
for high hx and small hi the motion in layer 1 is driven by the upper liquid. This is the case of regions (III) and 
(IV), where the direction of circulation close to the interface changes its sign. This is so because the buoyancy force 
acting on liquid 2 is so strong that it does not only drive the motion in the upper layer but it also modifies the flow 
direction near the interface. In (HI) two convection cells are found in liquid 1, but the one close to the lower surface 
is no longer driven by liquid 2. In region (IV) one observes the opposite scenario to (I), i.e. two convection cells but 
with the liquid flowing from the cold to the hot side near the interface. (The behavior of the lower liquid is analogous 
to that of a single liquid with 7 < HJ].) Contrary to the zero gravity problem states with one of the two layers 
at rest are not found in the present problem. 

Let us now examine the temperature field in the basic state (Fig. 3(b) I. The four regions do not coincide with 



these of the velocity field but they are related. Region (I) covers more than half the surface of Fig. 3(b) Liquid 2 
is unstably stratified (temperature increases with depth), in liquid 1 one finds a lower layer close to the wall which 
is unstably stratified, and above it one distinguishes a stably stratified sublayer (temperature decreases with depth). 
Region (II) is the only closed area: the lower liquid 1 is stably stratified whereas the upper liquid 2 is characterized 
by two stable stratified sublayers adjacent to the wall and the interface, in the middle there is an unstable stratified 
region, giving raise to a S-shaped profile. Region (III) shows the most complicated temperature profile: liquid 2 
exhibits two unstably stratified sublayers at the boundaries and a stable sublayer in the intermediate region, liquid 1 
consists in three sublayers, the intermediate being stably stratified while the two others are unstably stratified. Zone 
(IV) covers the low hi and almost all hx values, the global profile is similar to that of region (I), with the interface 
located at smaller hi so that liquid 1 is unstably stratified. 



IV. EVOLUTION EQUATIONS FOR THE PERTURBATIONS. 

As soon as a horizontal temperature difference AT is applied, convective cells typical of the basic state set in, 
but a further increase of AT may destabilize this basic fiow. To analyze the stability of the reference state with 
respect to infinitesimally small perturbations, let us write the general solution of the problem under the form u*-') = 
{uf + u'W,i;'W,m;'W), TW - T_ = x + + 6''W, = pW where (u'W, i;'^, w'W), 6l'W and p'(') denote 

the perturbations of velocity, temperature and pressure fields respectively. 

The perturbations are decomposed into a sum of normal modes (u',0',p') — (u'(z), 6''(z),p'(2;)) expi(k • x + wi), 
where oj denotes a complex frequency and k the wavenumber, with component in the streamwise direction and 
ky in the spanwise direction. The primes will be dropped for clarity in the ensuing equations. After eliminating the 
pressure in the momentum equation and omitting non linear terms, we are left with the following equations for the 
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perturbations, for the lower liquid, 1: 

9^4i(;(i) = (L(i) +/c2)a^,ii,(i) _L(i)/c2y;(i) +i?afc2g,(i) _j^^p^~i^(i)^^^^(i)^ (^^^^^ 

9,20(1) ^ (A:2 + z(cj + 4'^fc.))e^'^ +w^') +9^r(i)w(i), (11c) 
for the upper liquid, 2: 

= ^(2)^(2) + ^p^-la^u(2)^(2) ^ ^ / ^^^^^(2) „ ^(2) ^(2)\ ^ (i2a) 

1/9,4 ^ (i./fc2 + L(2))9,2 _ L<^^1 w^^^ + Ra k^ a 6^-^'^ - i k,Pr~'w<^^^ d,2U^o\ (12b) 

^9,20(2) ^ (Kfc2+i(c^ + 42)^^))0(2)^j2)_j_ 5^^(2)^(2)^ (^2C) 

wherein the following notation has been used 

L(i) = k^ + iPr-\uj + ul)^''k^), 
L(2) = uk^ + iPr-\u + u^Qh^). 

The corresponding boundary conditions are 



z = -1 ^ u(i) = u-^i) = 9,t(;(i) = 61(1) = 0; (13a) 

z = a^ u(2) = u;(2) = 9,17,(2) ^ ^(2) ^ Q. (13^) 

z = 0^ u(i)=77(2), 9,u;(i) =9,w;(2), u;(i) = u;(2) = 0, (13c) 
0(1) = 0(2), 9,0(1) = A9,0(2), /i9,w(2) - 9,u(i) = 7fc,Afa0(i), 
/792«;(2) -9>(i) = Afa/fc20(i). 

We are faced with an eigenvalue problem which is solved by means of a tau spectral method, by approximating the 
eigenf unctions with Chebyshev polynomials of order 16. The marginal curves are found by searching the set of values 
kx,ky, Ma for which the rate of temporal growth is zero. 

V. RESULTS OF THE LINEAR ANALYSIS 

Three mechanisms are able to destabilize the system under study. Basic temperature profiles in Fig. 3(b) present 



unstably stratified zones (temperature increases with depth) and other zones which are stably stratified (temperature 
decreases with depth), in the unstably stratified regions the Rayleigh-Benard instability may arise. A second mecha- 
nism is the usual Benard-Marangoni instability when the interface is colder than the adjacent fluid. Finally, a third 
destabilizing mechanism arises even when the core of the layers is colder than the interface, provided the velocity flow 
is strong enough to overcome the stabilizing vertical temperature field . 

Owing to the great number of parameters involved in a two-layer configuration, the results will be discussed in 
the case of selected pairs of liquids used in experiments. The calculations have been made for the pair 5 cSt silicon 
oil (upper liquid) and perfluorinated HT-70 (lower liquid), but later on we also consider the combinations: water 
(upper fluid) with perfluorinated hydrocarbon Fc-75 (lower liquid), and n-hexane (upper liquid) with acetonitrile 
(lower liquid). The parameter values of these fluids are gathered in Tab. U Dcformability is generally quantified by 

means of the so-called crispation number, defined as Cr = ^' ^J^ghW^ ' • shown in Tab. ^the values of Cr are very 
small for the pairs of liquids examined in the present work. Then the hypothesis of nondeformability of the interface is 
reasonable as it has been shown en earlier publications [2^l23|. Moreover, we will not pay attention to the particular 
hi values close to 0% or 100% because in these cases thin film motions (outside the scope of the present analysis) 
would appear. 
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A. Stationary longitudinal rolls 

First, we consider 2-D perturbations characterized hy — and to — (longitudinal stationary rolls). Our analysis 
extends for two layer systems that of Gershuni et al. |^ , who considered a one-single fluid layer bounded by horizontal 
rigid- free and perfectly thermal conductive surfaces, and Mercier and Normand "6], who considered the heat transfer 
at the free surface by means of a Biot number. 

Fig. 01 displays the critical Marangoni and the critical spanwise wavenumber versus the fractional bottom depth 
for several values of the total depth Ht (3,6 and 9 mm). For total depths of 3 and 6 mm, one observes three 
stationary branches, while for Ht — 9 mm an additional branch is clearly seen for low hi. Fig. 4(a) shows that 
the critical Marangoni number decreases with hx. This is due to the destabilizing effect of the increase of depth in 
Rayleigh-Benard instability. All the curves exhibit the most stab le zon e for hi « 60%. The jump from one branch to 
another is better appreciated by the discontinuities in ky in Fig. |4(b]1 The spanwise wavenumber ky for hx = 3 mm 



and hx — 6 mm increases with hi in the first branch, however it decreases in the second one together with that of 
hx — 9 mm, the dependence with hi in third branch is more involved for the three depths. In the branch jump the 
marginal Marangoni curve is bimodal 0, switching the absolute minimum between two wavenumbers. This modal 
change gives rise to different states in the two-layer system. 

In Fig. ISlare represented the isotherms and the velocity fields in the plane z — y as observed at the critical threshold. 
At small fractional bottom depths ( hi = 14% in Fig. |5(a)| ) strong temperature gradients are found in the bulk of 
liquid 1, whereas in liquid 2 they are concentrated near the interface z = 0, where they are less important. A similar 



behavior is exhibited by the velocity fields. For large fractional bottom depths ( hi = 80% in Fig. 5(b) I located at 
other branch the temperature gradients are concentrated in liquid 1 near the interface and do not penetrate into the 
bulk of the lower liquid, the velocity field is the largest in the lower liquid and is felt of roughly the same depth than 
the temperature gradients. Other stationary branches exhibit different kind of behaviors. 



B. Three dimensional perturbations 

However, one cannot exclude the possibility that the most unstable disturbances are three dimensional. If follows 
that a minimization process with respect to and ky is required to find the critical parameters. A typical stability 
surface is shown in Fig. |S1 In presence of lateral heating the reflection symmetry in the direction x is broken. Hence 
the marginal surface is only symmetrical under the changes ky — > —ky. The marginal surface is, in general, more 
complicated because several intermingled branches arise. The position of the minimum is given by the k^ and ky 
components of the critical wavenumber, which in the case of oscillatory instability also determine the direction of 
propagation ip of the waves. Recalling that the positive x axis is directed from the cold to hot side according to 
Eq. (O the hydrothermal waves will propagate from cold to hot for (p G [0°,90°], in contrast for Lp G [90°, 180°] the 
direction of propagation is the opposite. Note that the symmetry (y — —y) of the linear problem makes that the 
waves with wavenumber (kx,±ky) become unstable simultaneously. From now on, we will restrict to p d [0°, 180°] 
without lost of generality. 

In Fig. |7(a)| are represented the critical Marangoni number, frequency, angle and modulus of the wavenumber for 
the system HT-70 (lower liquid) and 5 cSt silicon oil (upper liquid) with hx = 9 mm. The critical Marangoni curve 
exhibits a local maximum A/1 and a local minimum at hi = 31%. The Marangoni number reaches its maximum value 
in M2, which is a crossing point between two branches. 

Beyond A/2, the behavior of the critical Marangoni is the same as in a one-layer system filled by the lower fluid. 
In such a case the critical Marangoni also decreases and the spanwise wavenumber increases with the depth Q . 

The position of M2 depends on the physical properties of both fluids. To study the effect of the transport coefficients, 
i.e. viscosity and thermal diffusivity k on the position of M2 we have considered an ideal system formed by two 
layers of 5 cSt silicon oil which differ only by their values of and k. In Tab. [H] are given the positions of A/2 for 
different values of the transport coefficients. By increasing the thermal diffusivity and viscosity of the lower liquid, 
the position of M2 is shifted towards higher values of hi; on the contrary if u and k are increased in the upper liquid 
the shift is towards lower values of hi. The location of M2 depends also on the geometry of the system: an increase 
of hx produces a shift towards a greater hi, however for high hx values the shift becomes saturated, as shown by the 
Marangoni curves of Figs. [3and|Hl 

Examination of the vertical profile of temperature of the basic state allow us to determine the positions of the local 
maximum Ml at 21% and of the local minimum at about 31%. Fig. 3(b) displays the vertical profile of temperature; 



it is seen that, for hx = 9 mm, starting from the region (I) and decreasing hi up to 31% the unstable vertical proflle 
of the upper fluid splits in two smaller unstable sublayers and a third one in the middle which is stable. The shape 
of the temperature proflle remains unchanged in liquid 1. 
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In the configuration represented by region (III), there is a greater extension of the stable regions than in (I), thus 
the scenario in (III) consists in an increase of Ma due to the progressive stabihzation of the upper layer when hi 
decreases. The border between regions (I) and (III) gives the position of the local minimum of Fig. |7(a)| However 
when region (IV) is attained the stratification in liquid 1 becomes unstable and the whole system is more and more 
unstable, with a decrease of Ma. This finds its roots in the two driving forces (thermocapillary and buoyancy), which 
are both destabilizing and cooperatives for the perturbations undergone by the lower liquid. As a consequence, the 
border between regions (III) and (IV) of Fig. |3(b]1 gives the position of Ml. Inside (IV) the decrease of Ma when hi 
is lowered reflects the greater importance of the destabilizing effect of the interfacial tension. The position of M 1 is 
predicted by the vertical temperature profile of the basic state and it will not be influenced by the Prandtl number. 

One distinguishes three main regions A, B, C in the curves of Figs. 13 and |H1 giving the critical parameters as a 
function of the relative depth of the liquids. 



1. region A 

The first region A extends from the lowest value of hi until M2. The critical Marangoni increases with the drop 
of the depth of the upper liquid because of the decrease of the buoyancy effects in it, indicating that the upper liquid 
dominates the dynamics in A. The pattern in this region consists of hydrothermal waves with zero or small angles 
of propagation. They propagate from cold to hot regions with an angle parallel to the gradient of temperature for 
/iT = 6, 9 mm, deviating to greater angles when the thermocapillary effects are relatively more important, as for 
/iT = 3 mm. Concerning the oscillatory modes, the angle of propagation is similar to single liquids with a high Pr 0. 
As seen in Fig. 7(d) region A is characterized by an increase of |k| with hi, with a small dependence of |k| on /ly. 



At /iT = 9 mm, where the buoyancy effects are more important, one finds a range of depth ratios 6% < hi < 16% 
where the critical modes take the form of stationary rolls. At the borders of this window, there is a codimension-two 
(cod-2) point formed by the interaction between a Hopf mode and a stationary mode of different wavenumbers. Such 
cod-2 points were also found in two-liquid layers heated from below p^ ]. 



2. region B 

Region B is characterized by a strong jump of the frequency at A/2 for /it = 6,9 mm, and weakly shifted by A/ii = 4 
for hx = i mm. For hx — 6,9 mm, M2 is characterized by a pair of Hopf modes with finite and different wavenumbers 
and frequencies which are simultaneously critical. At this point of cod-2 the configuration is the most stable with 
the greatest Ma. The difference between the wavenumbers of the two Hopf modes is due to the occurrence of a large 
spanwise component for h^ — 9 mm, while for hx — G mm there exists only a jump in the streamwise component. 
To our knowledge this kind of interaction between two oscillatory modes has not yet received attention in Benard- 
Marangoni problems. For /i^ = 3 mm, the presence of a small stationary branch near M2 prevents the existence of 
this kind of interaction, and the cod-2 point is generated by a stationary mode and a Hopf mode. Concerning the 
hydrothermal waves, the corresponding angles of propagation are greater than the angles of propagation in region A. 
The direction of propagation of the hydrothermal waves in this region is reversed, propagating from the hot to the 
cold side. The width of region B decreases with hx, from Ahi = 18% for hx = 9 mm to A/ii = 14% for /it = 3 mm. 
The modulus of the wavenumber varies in a more complicated way than in the region A, and is more sensitive to the 
total depth. In region B, the strong competition between the two layers is responsible for the oscillatory convection. 



3. region C 



The transition from the second region B to the third region C is smooth with a monotonous decrease of the frequency 
of the oscillatory branch. The third region is characterized by stationary longitudinal rolls and extends up to the 
highest hi. It spreads towards lower hi when the total depth is increased. The dynamics of this region is dominated 
by the lower layer. The qualitative behavior of the modulus k of the wavenumber is practically unaffected by the 
overall depth. 



In Fig. 



9(a) are shown the critical parameters for the bilayer formed by perfluorinated hydrocarbon Fc-75 (lower 
liquid) and water (upper fluid), whereas in Fig. 9(b) are reproduced the results for the couple acctonitrile (lower 
liquid) and n-hexane (upper liquid). Both setups have been experimentally studied in the case of a vertical heating 
[Hill El- The instability thresholds are evaluated for a total depth of 6 mm in both configurations. The critical 
Marangoni curves exhibit maxima at Ml and M2 as discussed above. In the case Fc-75/water, hydrothermal waves 
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which propagate from cold to hot in the region located between Ml and M2 are found, for the other values of the 
fractional bottom depth stationary longitudinal rolls are forecast. The system acetonitrile/n-hexane presents the same 
three regions as exhibited by the pair 5 cSt silicon oil/HT-70 with a cod-2 point created by two oscillatory modes 
with different frequency and wavenumbers. The general features discussed earlier in this section remain valid for all 
the cases which have been examined. 

Tuning the relative depths of both layers, according to our theoretical study it would be possible to obtain in the 
same experimental configuration with a fixed pair of liquids all the patterns observed in different single layer systems: 
from longitudinal stationary rolls to hydrothermal waves with directions of propagation ranging from small angles of 
propagation (as in one- liquid systems with high Pr) to high angles (as in one- liquid systems with small Pr). 

VI. CONCLUSIONS 

We have investigated the onset of convection in systems constituted by horizontal liquid layers subjected to a lateral 
heating. The system is bounded by two rigid and thermally conductive horizontal walls and the interface is supposed 
to be nondeformable. The acting forces are the interfacial tension and the gravity. They are cooperative in the lower 
layer and competitive in the upper one. The number of parameters is greatly enhanced with respect to the case of 
one-layer problems. We have considered three different combinations of liquids, namely bilayers of perfluorinate HT-70 
and 5 cSt silicon oil, perfluorinate Fc-75 and water, and acetonitrile and n-hexane; these fluids have been considered 
because they have already been the subject of experiments involving a vertical heating [l5l |. It is our hope that the 
results of the present paper will motivate researchers to repeat these experiments with a lateral heating. 

As soon as a horizontal temperature gradient is applied, a global circulation takes place in each layer. The tem- 
perature and velocity fields of the basic state have been derived and discussed. For a given experimental setup and a 
fixed heating, the two control parameters are the total depth and the depth ratio between the liquids. In a diagram 
representing the depth ratio versus the total depth, four regions for the temperature and velocity profiles have been 
distinguished. The main information drawn from this diagram can be summarized as follows. The interaction between 
interfacial tension and buoyancy can give rise to: 1) a counter-rotating cell in each layer, 2) two counter-rotating cells 
in the upper liquid and one cell in the lower one, in both cases one observes a flow from hot to cold near the interface ; 
3) to two counter-rotating cells in the lower liquid and a cell in the upper one, 4) a counter-rotating cell in each liquid, 
but with the flow being dragged near the interface from cold to hot in the last two cases. We have also determined 
the different temperature proflles that result from the interplay between thermal diffusivity and fluid motion. 

The linearized evolution equations for the perturbations of the basic state have been established and the ensuing 
eigenvalue problem has been solved. Bidimensional and three dimensional perturbations have been respectively 
considered. A physical interpretation of the different regions found has been proposed. The interest of the present 
analysis with respect a one-layer system is that it exhibits a great variety of behaviors depending on the nature of the 
liquids the total depth and relative depth between the layers. Among the most interesting results can be mentioned 
the occurrence of a cod-2 point arising as a consequence of the interaction between two Hopf modes with different 
frequency and wavenumber. Furthermore, three different patterns are predicted: 1) hydrothermal waves propagating 
from cold to hot side with a small angle, which is typical of lateral heating in one-liquid systems with high Pr, 2) 
hydrothermal waves propagating from hot to cold, and 3) longitudinal stationary rolls. 

The variety of results confirms that the problem of lateral heating in two- layer systems is a promising and interesting 
area of research, both from the theoretical and experimental points of view. It would be specially interesting to study 
the system close to Af2 when we have interplay between two oscillatory modes. 
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APPENDIX 

Expressions of the velocity and temperature profiles of the basic state in dimensionless units: 
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u, 



- ^ [-12aA/at/(l + 4z + 3z2) (14) 



,(2) 



48 (a + ^) f 

+ Ra {a'^ a fj, — ai' + 4 fi (a"^ a + i') z + 3 [a^ a fi + 3 av + A fiiy) z'^ + 8 {a + fx) z'^}] 
1 



[-12Maiy{a'^ -Aaz + 3z'^) (15) 

(16) 



iSa {a + ^) v 

+ Ra afi^a^iy + Aa (a'' a + z — 3 (4 a"^ a + 3 a + i^) + 8 a a {a + ^) z^}] 
1 



2880 (a + A) (a + /x) z/K ' 
60 Ma \^c? V (a A + k) + a A (a^ — k) z — 6 a (a + A) v nz^ 



\a (a + A) ;/ K z"^ — 3 a (a + A) // k z*} 



+i?a |at^ (5a'^A + 9aK + 4/iK) — a'* a (4a^A + 9aA/i + 5/iK) 

- (i^ (-5 + 9 a K + 4 /i k) + a (4 a'' + 9 ^ - 5 ^ k)) A z 

+30 a (a + A) [c? 0^1- v) k z^ + 40 (a + A) /i (a^ a + i') k z^ 

+ 15 (a + A) (a^ a /I + 3 a i/ + 4 ^ z^) k z'' + 24 (a + A) (a + ^) v k. z^ }] 

r(2) = — [ (17) 

2880 a (a + A) (a + /i) i/K ^ ^ ' 

60 Ma {- a^ (a A + k) + a^ (-a^ + k) z + 6 a^ (a + A) j/z^ - 8a (a + A) z^ 

+3 (a + A) j/z"*} +i?a {a^ (- (z^ (5 a^ A + 9 a k + 4 ^ k)) + a^ a (4 A + 9 a A /i + 5 ^ k)) 

+a (i/ (-5 + 9 a K + 4 k) + a (4 a*" + 9 a^ ^ - 5 a^ k)) z - 30 a^ (a + A) {(? a[i~v) z^ 

-40 a (a + A) {c? a + j/) z^ + 15 (a + A) {(? a (4 a + 3 ^) + ;/) z'' - 24 a a (a + A) (a + ^) z^ }] 

Note that the lower fluid (superscript 1) is defined in the interval z e [—1,0] and the upper fluid (superscript 2) in 
the interval z S [0, a], according to the definitions given in section Ull 
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FIG. 1: Diagram of the two-layer system under study. Two superposed immiscible liquids bounded horizontally by conductive 
rigid walls are subject to a horizontal gradient of temperature. The interface is supposed to be non deformable. The gravity 
and interfacial tension are the forces acting on the system. 




FIG. 2: Driving forces: gravity (left) and interfacial tension (right). The buoyancy and thermocapillary forces act in the same 
direction in the lower layer, and in opposite direction in the upper layer where they are competing. 



14 




15 




FIG. 4: Critical Marangoni (left) and critical wavenumbcr (right) for 2-D longitudinal perturbations {kx = 0) as a function of 
the bottom depth ratio, for three different values of the total depth. The results correspond to the pair HT-70 (lower liquid) 
and 5 cSt silicon oil (upper liquid). 
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FIG. 5: Isotherms (above) and velocity fields (below) in the plane z-y are represented at critical thresholds for 5 cSt silicon 
oil/HT-70 system and hr = 9 mm, for two bottom depth ratios corresponding to different branches: hi = 14% (right ), 
hi =80% (left). 




FIG. 6: Surface of neutral stability for three dimensional perturbations. In this example hr = 6 mm, and hi = 78%. The 
breaking of kx, ky symmetry is clearly appreciated. The stability threshold is located at k = (4.2, 7). 
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FIG. 7: Critical paxameters for the system 5 cSt silicon oil (upper liquid) and HT-70 (lower liquid) with total depth 9 mm. (a) 
Critical Marangoni, (b) critical frequency, (c) angle of the critical wavevector and (d) modulus of the critical wavevector. Solid 
line: oscillatory behavior, dashed line: stationary behavior. In the points where there is an interaction between two branches 
are displayed their extensions in the Marangoni curve. Different characteristic patterns are obtained in regions A, B and C ( 
see the text). 




FIG. 8: Critical parameters for the system HT-70 (lower liquid) and 5 cSt silicon oil (upper liquid) with total depth 6 
and 3 mm (b). Solid line: oscillatory behavior, dashed line: stationary behavior. 
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FIG. 9: Critical parameters for the pair Fc-75 (lower liquid) and water (upper liquid) (a) and for the pair acetonitrile (lower 
liquid) and n-hexane (upper liquid) (b). The total depth for both configurations is /it = 6 mm. 



21 





P 


V A 


Cp 


a 


7 


Jr^T 
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[kg/m ) 




rn^/s) {J/ms 


K) {J/kgK) (10-3 i^-i) 


(JN/mKj 






Silicon oil 5 cSt 


920 


5 


0.117 


1590 


1.05 




62.512 




HT-70 


1680 


0.5 


0.07 


962 


1.1 




11.54 




Silicon oil 5 cSt /HT-70 


0.548 


10 


1.671 


1.653 


0.954 


-7.3 -10-5 




2.1 -10-® 


Water 


997 


0.893 


0.609 


4180 


0.257 




6.111 




c C~ i 


1 7«C\ 
1 ( ou 


0.945 


0.063 


1046 


1.4 








Water/Fc-75 


0.566 


0.945 


9.59 


3.996 


0.183 


-4.7-10"^ 




3.3-10-'' 


n-Hexane 


655 


0.458 


0.12 


2270 


1.41 




5.675 




Acetonitrile 


776 


0.476 


0.118 


2230 


1.41 




4.381 




n-Hexane /Acetonitrile 


0.844 


0.962 


1.017 


1.018 


1 


-1 -10-^ 




7.8-10-'' 



TABLE I: Parameter values of the three pairs of liquids studied. The parameter ratio in each different configuration is also 
given. The Crispation number Cr is given for a depth of /i'^^ = 3 mm. 
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V 








0.1 


84% 


0.1 


90% 


1 


81% 


1 


81% 


10 


68% 


10 


62% 



TABLE II: Position of M2 as a function of the viscosity ratio (loft) and thermal viscosity ratio (right) of the two hquids. The 
case V = 1, K = 1 corresponds to an ideal configuration formed by two layers of the same fluid (5 cSt silicon oil). 



